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ABSTRACT 

Observational evidences point to a common explosion mechanism of Type la supernovae based 
on a delayed detonation of a white dwarf. Although several scenarios have been proposed and 
explored by means of one, two, and three-dimensional simulations, the key point still is the 
understanding of the conditions under which a stable detonation can form in a destabilized white 
dwarf. One of the possibilities that have been invoked is that an inefficient deflagration leads to 
the pulsation of a Chandrasekhar-mass white dwarf, followed by formation of an accretion shock 
around a carbon-oxygen rich core. The accretion shock confines the core and transforms kinetic 
energy from the collapsing halo into thermal energy of the core, until an inward moving detonation 
is formed. This chain of events has been termed Pulsating Reverse Detonation (PRD). In this 
work we explore the robustness of the detonation ignition for different PRD models characterized 
by the amount of mass burned during the deflagration phase, Mdefi- The evolution of the white 
dwarf up to the formation of the accretion shock has been followed with a three-dimensional 
hydrodynamical code with nuclear reactions turned off. We found that detonation conditions are 
achieved for a wide range of M^ c ^. However, if the nuclear energy released during the deflagration 
phase is close to the white dwarf binding energy (~ 0.46 x 10 51 erg =>■ M^ofl ~ 0.30 M Q ) the 
accretion shock cannot heat and confine efficiently the core and detonation conditions are not 
robustly achieved. 

Subject headings: Supernovae: general - hydrodynamics - nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

Type la supernovae (SNIa) are the most en- 
ergetic transient phenomena in the Universe dis- 
playing most of its energy in the optical band 
of the electromagnetic spectrum, where they ri- 
val in brightness with their host galaxies during 
several weeks. The importance of SNIa in as- 
trophysics and cosmology is highlighted by their 
use as standard (or, better, calibrable) candles to 
measure cos mic distances and related cosmological 
parameters dPerlmutter et al . 19981 ISchmid~ et al. 
Il998i : iRiess et all l200lh . However, in order to 
achieve a high precision in the distance deter- 
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mination as required, for example, to determine 
the equation of state of the dark energy com- 
ponent of our Universe, it is necessary to un- 
derstand the physics of SNIa explosions. From 
the theoretical point of view, the accepted model 
of SNIa consists of a carbon-oxygen white dwarf 
(WD) near the Chandrasekhar mass that accretes 
matter from a companion in a close binary sys- 
tem. This model accounts for the SNIa sam- 
ple homogeneity, the lack of prominent hydrogen 
lines in their spectra, and its detection in ellip- 
tical galaxies. Massive WDs are extremely un- 
stable bodies in which a modest release of en- 
ergy can produce a huge expansion (i.e. their 
ratio of binding to gravitational energy is only 
~ 15%, while in a normal star it is ~ 50%). There 
are two main ingredients of the standard model 
that are still poorly known: the precise configura- 
tion of the stellar binary system and its evolution 
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prior t o thermal runaway o f the WD (lLanger et al 



2000; Nomoto ct al. 2003; Piersanti et al. 2003; 



Han fc Podsiadlowskil 12004 iBadenes et alJ 12007 : 
Hachisu et al. Il2008h. and the ex plosion mechanism 
(Hille brandt fe Niemeverl l2000h . Fortunately, as 
long as a carbon-oxygen WD reaches the Chan- 
drashekar mass, its previous evolution is not crit- 
ical for the explosion because the WD structure 
is determined by the state of a degenerate gas of 
fermions that can be described with a single pa- 
rameter, the mass of the star. This fact leaves 
the explosion mechanism as the most relevant un- 
known concerning SNIa. 

In spite of continued theoretical efforts dedi- 
cated to understand the mechanism behind SNIa, 
realistic simulations are still unable to provide a 
satisfactory description of the details of these ther- 
monuclear explosions. Nowadays, there is a con- 
sensus that the initial phases of the explosion in- 
volve a subsonic thermonuclear flame (deflagra- 
tion), whose propagation competes with the ex- 
pansion of the WD. After a while the corrugation 
of the flame front induced by hydrodynamic in- 
stabilities culminates in an acceleration of the ef- 
fective combustion rate. However, the huge dif- 
ference of lengthscales between the white dwarf 
(~ 10 8 cm) and the flame width (< 1 cm) prevents 
that large-scale numerical simulations resolve the 
deflagration front, making it necessary to imple- 
ment a subgrid model of the subsonic flame. The 
precise value of the effective combustion rate is 
currently under debate, as it depends on details of 
the flame model. Recent three-dimensional calcu- 
lations by different groups have shown that pure 
deflagration models always give final kinetic en- 
ergies that fall short of 10 51 ergs, while leaving 
too mu ch unburnt carbon and oxygen close to the 
center jReinecke et al l 120021 IGamezo et alJ 120031: 
Garcia-Senz fc Bravol l2005t iRopke et alJ l2007af ). 
Because both signat ures are at odd s with obser- 
vatio nal constraints, IGamezo et al. (120031 I2004L 
20051 ) concluded that the only way to reconcile 
three-dimensional simulations with observations is 
to assume that a detonation ignites after a few 
tenths of a solar mass have been incinerated sub- 
sonically (delayed detonation). 

The idea that a delayed detonation is some- 
how involved in thermonuclear explosions lead- 
ing to SNIa h as be en around for many years. 
Ivanova et al. ( 1974 ) proposed that a delayed 



detonation might take place after a long-range 
pulsation of the white dwarf, in what is known 
as the Pul sating D e layed Detonation (PDD) 



scenario. 



Khokhlovl (|l99lh carefully computed 



one-dimensional delayed detonation models that 
were successfully compared with observations of 
SNIa explosions ( Hofhch et al.l 19951: Howell et al 



Fesen et al 
20061 



2007 1) and th e ir rem nants (jBadenes et al 
In these one- 



2008; iRest et al.l [2008) 



20011: iQuimbv et all 120071: iGerardv et alJ 120071: 



dimensional models the detonation initiation had 
to be postulated, as there was not identified any 
sound physical mechanism by which such a tran- 
sition could hap pen in an unconfined white dwarf 
(|Niemeverlll999T) . 

Ignition of a detonation in a white dwarf 
can happen in two basic ways, either through 
a deflagration-to-detonation transition (DDT) or 
because of a sudden energy release in a confined 
fluid volume (hereafter: confined detonation ig- 
nition or CDI). The essential feature of detona- 
tion initiation is the formation of a non-uniformly 
preheated region with a level of fluctuations of 
temperature, density, and chemical composition 
such that a sufficiently large mass burns be- 



fore a sonic wave can cross it (Khokhlov 1991 



Niemever fc Wooslevlll997h. The therma l gradient 
needed is dBlinnikov fc Khokhlov 1987 : Wooslev 
19901 : lKhokhlovlll99lh : 



VT < 



sound T\ 



(i) 



where A is a numerical coefficient, t\ = T/T is the 
induction timescale at temperature T, and ~ 
0.04 — 0.05 is the Frank-Kamcncctskii factor: 



e = - 



<91nT 
d In 71 



(2) 



Such fluctuations could be produced by a vari- 
ety of mechanisms: adiabatic pre-compression in 
front of a deflagration wave , shock heating, mixing 
of hot ashes with fresh fuel |Khokhlovlll99ll) . accu- 
mulation of pressure waves due to a topologically 
complex geometrica l struc ture of the flame front 
(|Wooslev fc Weaverl I1994T ). or transition to the 
distri buted burning regime (jNiemever fc Wooslev 



1997). Among the proposed mechanisms of DDT, 



turbulence pre-conditionin g has r eceived the most 
attention. iKhokhlov et al.l (|l997l ) determined cri- 
teria for a DDT in unconfined conditions, such 
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as those realized during the expansion of a white 
dwarf following a deflagration phase. For a DDT 
to be feasible, the turbulent velocity has to ex- 
ceed the laminar flame velocity by a factor < 8 
at a lengthscale comparable to the detonation 
wave thickness. This criteria was fulfilled for 
flame densities in the range 5 x 10 6 g/cm 3 < 
p < 2 — 5 x 10 7 g/cm 3 for reasonable assump- 
tions. At densities in excess of 10 s g/cm 3 a DDT 



is quite unlikely ([Khokhlov et al.1 119971 but see 



Zingale fe Dursil who pointed to bubbles fragmen- 



tation as a way to increase the flame surface and 
facilitate a DDT at p ~ 2 x 10 8 g/cm 3 ). 

Small-scale simulations are needed to ascertain 
if the necessary conditions for a DDT are ac- 
tually achieved during a white dwarf explosion 
driven by a deflagration wave. Up to the present, 
such studies seem to disfavor a DDT in view of 
the robustness of subsoni c flames aga i nst in ter- 



action with vortical flow ( Ropke et al. 2004 ) for 



the maximum expe cted velocities at the integral 
scale, ~ 100 km/s dLisewski et al.l [2000h . More- 



Niemeyer (1999) estimated that the size of 



over, 

the fluctuations that can be expected from turbu- 
lence at the critical densities is ~ 3 orders of mag- 
nitude smaller than required for a DDT to happen. 

Three-dimensional numerical simulations of 
white dwarfs undergoing a slow defla gration have 



provid ed new ways to achieve a CDI. iPlewa et al 



(|2004f ) followed the evolution of a white dwarf as- 
suming ignition in a single slightly off-centered 
bubble, and found that a CDI might result from 
the convergent flow at the antipodes of the ini- 
tial ignition point. In their model, the det- 
onation was possible because of the confining 
gravitational field, hence it was termed Gravi- 
tationally Confined Detonation (GCD, see also 



Ropk e et al 
20071 



200714 
et all I 



Plewal 120071: iTownslev et al 
Jordan et al.1 120081 ). iBravo fe Garcia-Senz 
(|2005 ) described another CDI scenario in which a 



detonation was triggered by inertial confinement 
due to an accretion shock born around a carbon- 
oxygen white dwarf core after pulsation caused by 
burning of ~ 0.1 — 0.3 M Q : the Pulsating Reverse 
Detonation (PRD0- In this work, we explore the 
robustness of detonation ignition conditions in the 



A conceptually similar m odel was proposed 
by iDu nina-Barkovskava et a k| d200lT) . See also 



PRD scenario, as a function of the mass burned 
during the deflagration phase prior to pulsation, 
Afdcfi- The evolution of the white dwarf up to the 
formation of the accretion shock has been followed 
with a three-dimensional hydrodynamic code with 
nuclear reactions turned off in order to determine 
the most extreme conditions achieved by the white 
dwarf during the re-collapsing phase of the pulsa- 
tion. In the next section we summarize the main 
features of the proposed scenario, and compare 
PRD models with one-dimensional PDD models. 
After that, we analyze if detonation conditions 
are achieved in our PRD models. We then discuss 
the implications that the presence of a small cap 
of helium on the white dwarf at the moment of 
thermal runaway might have for CDI. Finally, we 
present the conclusions of this work. A companion 
paper (Bravo et. al., 2009; hereafter paper II) is 
devoted to present the detailed evolution of PRD 
models after detonation ignition and their final 
outcomes. 

2. Briefing of the pulsating reverse deto- 
nation model 

In the PR D scenario the explosion pro ceeds in 
three steps ( Bravo fc Garcia-Senz 20061) : 1) an 
initial pre-conditioning phase of the WD whose 
result is an expanded structure with a positive 
gradient of the mean chemical weight, 2) forma- 
tion of an accretion shock that confines the fuel 
volume and, 3) launch of a reverse, converging, 
detonation wave. We have probed the PRD con- 
cept by performing three-dimensional simulations 
of the explosion. The initial model consists of a 
1.38 Mq WD made of carbon and oxygen. The 
hydrodynamic evolution starts with the ignition 
of a few sparks and burns < 0.3 M©, mostly to 
Fe-group elements, during the first second, releas- 
ing < 5 x 10 50 ergs of nuclear energy. The released 
energy resides for the most part in the outermost 
layers leading to WD pulsation. 

The details of the first phase, that spans 
the first second after thermal runaway, have 



been known for some time dPlcwa et al. 2004; 



Dunina-Barkovskaya & Imshcnnik (2003) 



Garcia-Senz fc Bravol 20051 : Livne et al. 2005 ) so 
we will only briefly sketch here the evolution dur- 
ing this phase. Even though the precise configura- 
tion at thermal run away is difficult to determine 
(jHoflich et al.ll2003l h current works suggest a mul- 
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tipoint ignition in which the first sparks are lo- 
cated slightly off-cen t er (IGarcia-Senz fc Wooslev 
1995[ IWooslev et"aT1 120041 : IWunsch fc Wooslev 
2004). If the number of sparks is too small or 
the burning rate is too slow the nuclear energy 
released is not eno ugh to unbind the star and 
the explosion fails (jRopke et al.ll2007b). This is 



know n as the "bubble catastrophe" (jLivne et al 
2005). The bubbles float to the surface before 



the combusti o n wave can propagate substantially 
(jPlewa et all 12004 IGarcia-Senz fc Bravo! 120051 ) 
and the star remains energetically bound. This 
behavior produces a composition inversion, i.e. 
the internal volume is plenty of fuel, cold carbon 
and oxygen, while the ashes of the initial com- 
bustion, mostly hot iron and nickel, are scattered 
around. 

The second phase of the explosion starts when 
the deflagration quenches due to expansion and 
ends when an accretion shock is formed (several 
seconds after initial thermal runaway) by the im- 
pact of the in-falling material onto the carbon- 
oxygen core. An important feature is that the 
decompression and further expansion of the cold 
core in his way to regain hydrostatic equilibrium 
imparts mechanical work to the iron-rich atmo- 
sphere. The core-atmosphere energy transfer con- 
tributes to the strength of the shock formed when 
the atmosphere re-collapses. The evolution from 
the time of maximum expansion up to the mo- 
ment of formation of the accretion shock can be 
compared to simulations of the self-gravitating col- 
lapse of a cold gas sphere, a well-known hydrody- 
namica l test that can be found, for instance in 
Evrardl (1988). The structures seen in the third 
column of their Fig. 5 qu alitatively agree with 
those appearing in Fig. 1 of iBravo fc Garcia- Senz 
(2006) in spite of the different initial conditions. 



2.1. Numerical simulations of the pulsat- 
ing phase 

We have simulated the evolution of the white 
dwarf during the second phase of the PRD explo- 
sion mechanism with a Smoothed Particle Hydro- 
dynamics (SPH) code who se main features have 
been described elsewhere (|Garcfa-Senz fc Bravo 
2005). The number of particles used in the exper- 
iments was N = 250 000, which provided a good 
spatial resolution, as detailed later. The initial 
model for these numerical experiments is given by 



the output obtained at the end of the deflagra- 
tion p hase computed as in Garcia- Senz fc Bravol 
(|2005l ). Our working hypotheses are that dur- 
ing the deflagration phase a mass in the range 
Mdefl ^5 0-30 M Q is burnt and that the final out- 
come of the explosion is determined by the value 
of Mdefl- We disregard any effect due to differ- 
ent possible deflagration histories that might lead 
to a given Mdefl with different chemical composi- 
tions (e.g. different proportions of Fe-group and 
intermediate-mass nuclei) and chemical profiles. 
Hereafter, the present models will be designated 
with the acronym DFnn, where 'nn' stands for the 
hundredths of solar mass burned during the de- 
flagration phase, for instance model DF11 means 
that Mdefl = 0.11 M Q . 

In Fig. Q] there is shown the evolution of models 
DF11, DF18 and DF29 (characterized by M dofl = 
0.11 M , M de a = 0.18 M Q , and M de fi = 0.29 M , 
respectively) through the deflagration and the pul- 
sating phases. These models started from differ- 
ent number of igniting bubbles and used differ- 
ent flame velocities with the result that a range of 
M defl = 0.11 - 0.26 M Q was obtained. The pro- 
cess of generati ng the initial model w a s the same 
as described in Garcia-Senz fc Bravo j2005l ). the 
location of the center of mass of the bubbles and 
their incinerated mass are given in Table [T] Mod- 
els DF11 and DF18 used only the first six bubbles 
in the Table, while model DF29 started from the 
seven bubbles listed. The flame velocities were: 
Udcfi = 100, 150, and 200 km s" 1 for models DF11, 
DF18, and DF29 respectively. 

As can be seen in Fig. [TJ most of the thermonu- 
clear burning took place during the first second for 
all models. Two factors determined the quench- 
ing of the burning (in this context by quenching 
we mean that the nuclear timescale became much 
larger than the dynamical timescale): global den- 
sity decrease due to white dwarf expansion, and 
bubbles migration to the external layers. The sim- 
ulations described hereafter belong to the evolu- 
tion following maximum expansion of the white 
dwarf. From that time on the evolution was fol- 
lowed with the nuclear reactions switched off in 
order to expose the fuel to extreme conditions dur- 
ing the collapse and check if a stable detonation 
could be achieved. Note that this is nothing but 
a convenient approximation to the real situation 
in which some preliminary burning is expected 
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before maximum compression, thereby modifying 
the maximum density and temperature achieved 
during pulsation, as shown in paper II. The three- 
dimensional evolution of the bubbles during the 
deflagration phase of model DF18 is shown in 
Fig. [2] The first noticing feature is that the evolu- 
tion of the bubbles is similar to each other: all of 
them grow in mass through burning and rise to the 
surface of the white dwarf at the same time that 
their sha pe becomes toroidal due to shear forces o n 
their top (jZingale et alj|2005l : lllopke et al.ll2007bl ). 
At t = 0.66 s fuel takes up the whole central vol- 
ume of the white dwarf. At t = 1.05 s the flame 
has virtually quenched and the bubbles have ex- 
panded laterally to completely surround the core 
of the star. Afterwards, the ashes hide the evolu- 
tion of the inner fuel during the pulsation. 

The evolution of model DF18 during the first 
pulsation, after deflagration quenching, is shown 
in Figs. [3] and HJ and the structure of the white 
dwarf at the time of accretion shock formation, 
t acc , is shown in Fig. [5l The inner ~ 1.10 M© 
pulsate in phase and reach a state of maximum 
compression at t = 7.16 s, at which time they 
make up a compact nearly hydrostatic core. At the 
same time, the outermost layers continue expand- 
ing with velocities larger than the escape velocity, 
while the intermediate layers continue falling onto 
the core with velocities as large as 5 000 km s _1 . 
The successive stretching and shrinking of the C- 
O rich core can be better seen in Fig. 0J Although 
the chemical structure carry the imprint of the 
bubbles evolution and lacks spherical symmetry, 
but for the central volume, the fuel mass fraction 
at radius below ~ 5 000 km is high: Xf ue i > 0.4. 

The mechanical and thermal structure at t acc 
display a high degree of spherical symmetry, as 
shown by Fig. [5] The picture of the radial velocity 
reveals a large gradient above r ~ 3 000 km. The 
opposed gradients of radial velocity and density 
lead to a quite homogeneous impact pressure (due 
to matter infall onto the core) as high as P ram = 
pv 2 ~ 10 23 dyn cm -2 in the range of distances r ~ 

3 000 — 5 000 km. This impact pressure takes over 
the thermal pressure above r ~ 4 000 km, so that 
the largest deceleration and subsequent heating of 
the in- falling matter takes place in a ring (in fact, 
a narrow shell in three dimensions) at r ~ 3 000 — 

4 000 km, as seen in the temperature plot. The 
fuel mass fraction in these locations at i acc lies in 



the range X\ 



fuel 



0.4-0.8. 



2.2. Comparison between PRD and PDD 
models 

In order to fully understand the implications of 
the three-dimensional calculations of the pulsating 
phase of the supernova it is convenient to compare 
the results of the present PRD simulations with 
those obtained with a one-dimensional code (PDD 
models). The key difference between the simula- 
tions performed in one and three dimensions is 
the impossibility of the former to change the se- 
quence of the mass elements. This has two im- 
portant consequences. First, in one-dimensional 
calculations the material burned close to the cen- 
ter never catches the fuel located ahead of it. 
Second, in one-dimensional calculations the un- 
burned mass always lies above the flame blocking 
the expansion of the underlying layers and facil- 
itating a more efficient burning. In contrast, in 
three-dimensional models the floatability of ignit- 
ing bubbles allows them to reach the external lay- 
ers with only a moderate transfer of mechanical 
work to their surroundings. Thus, in three dimen- 
sional PRD calculations the ashes retain most of 
their specific energy (kinetic and thermal), which 
enables the formation of a robust accretion shock 
that ultimately triggers a detonation. In contrast, 
in the PDD mechanism, the detonation is thought 
to be due to compression of the mixing layer be- 
tween unprocessed and processed material, which 
is believed to grow only if t here is a long-rang e 
pulsation of the white dwarf ( Hoflich et al.lll995h . 

In Figs. [5] and [7] there are represented the 
angle-averaged specific energy and chemical pro- 
files, respectively, of the three-dimensional model 
DF18 (top panel) and a one-dimensional calcula- 
tion of a pulsating delayed d etonation model (bot - 
tom panel: model PDDe in lBadenes et al.l l2003). 
In the three-dimensional calculations the energy 
resides for the most part in the outer 0.15 M Q . 
Hence, the pulsating motion of this material is de- 
coupled from the rest of the structure. In contrast, 
in the one-dimensional calculation the energy is 
concentrated in the inner 0.2 M©, that behaves 
like a piston driving the pulsation of most of the 
structure nearly in phase. In the PDD mecha- 
nism the detonation is initiated at the edge of the 
incinerated core during the collapsing phase and 
propagates outwards through the in-falling matter 
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(|Khokhlov et al.lll993l ). 

An important requirement of the PDD mecha- 
nism is the necessity to burn a quite precise mass 
of fuel, ~ 0.3 M Q , during the initial deflagration 
phase in order to achieve a large amplitude pul- 
sation (necessary to create a wide enough mixed 
layer) yet not unbind the white dwarf. This is not 
the case in the PRD models, because: 



evolution of the velocity profiles of the PDD model 
is displayed. Except for a few layers at the top of 
the white dwarf, the star pulsates in phase without 
development of any shock wave (unless a detona- 
tion is assumed to happen). These calculations 
show that the chemical composition inversion is 
a crucial ingredient in the evolution of the white 
dwarf up to the moment of detonation ignition. 



1. The detonation is not launched as a result 
of turbulent mixing of fuel and ashes, thus 
avoiding the necessity to achieve a large am- 
plitude of the pulsating motion. 

2. The relative amplitude of the pulsation is 
not uniform for all the mass layers of the 
white dwarf. The external layers experience 
a large amplitude because they are rich in 
hot ashes, freshly synthesized in the float- 
ing burning bubbles, while the 1 Mq core 
remains at moderately high density during 
the pulsation. Thus, in the PRD scenario a 
large amplitude of the external layers is ob- 
tained irrespectively of the precise amount 
of mass burnt during the deflagration phase. 

Finally, we compare the evolution during the 
pulsation of the structures obtained at the end of 
the deflagration phase for the three-dimensional 
model DF18 and that of a one-dimensional model 
constrained to burn subsonically the same mass as 
DF18. First, we have mapped the DF18 model 
at 4.2 s into an angle-averaged version and its 
posterior evolution has been followed with a one- 
dimensional hydrocode0 with the nuclear reactions 
turned off. A sequence of the velocity and entropy 
profiles obtained this way is shown in Fig. [5] In 
this Fig. it can be appreciated clearly the process 
of formation of the accretion shock, whose main 
distinctive feature is the sudden increase of the 
specific entropy. Second, we have computed a one- 
dimensional PDD model with Md e fl = 0.18 M Q . 
By construction, the difference between this PDD 
model and the one-dimensional version of DF18 re- 
sides in the distribution of fuel and ashes, as well 
as the specific energy they transport. The con- 
sequences are plain to see in Fig. |H1 in which the 



2 The one-dimensional hydroco de, based on the finite dif- 
ferences scheme pr oposed by | Co lgate &:~Whitc ( 19 63), is 
the on e described in lBadenes et alTl2003ft" and lBravo et aT] 
(1993ft 



3. Confined detonation ignition 

We now turn to the question of the robustness 
of the CDI in the PRD scenario. First we examine 
the physical conditions in the accretion shock, as 
obtained in three-dimensional simulations of the 
pulsating phase. By comparing these conditions 
with suitable detonation criteria, we then analyze 
the most favorable conditions for detonation igni- 
tion in our three-dimensional simulations. 

3.1. Physical conditions of shocked matter 

The location and physical conditions at the ac- 
cretion shock are very sensitive to the amount of 
mass burned during the initial deflagration phase. 
A summary of the properties of the accretion 
shock for three calculated models is given in Ta- 
ble [2j The mass at which the accretion shock 
is formed, M acc , is larger for smaller Mdofl be- 
cause less energy is available to power the ex- 
pansion of the white dwarf. The relationship be- 
tween both masses can be approximated in the 
range of M defl explored by: M acc w -4.18Mj cfl - 
1.22Mdofl + 1.50, where all masses are expressed 
in M Q . Other quantities that can be fitted to an 
analytic function of Af acc are: the time of forma- 
tion of the accretion shock, i acc w — 15.9A/ acc + 
25.1 s, and the central density at that time, p c $ rs 
exp (8.87M£. C - ll.OMacc + 1.50), where p c8 is in 
units of 10 s g cm~ 3 . 

Once formed, the accretion shock remains con- 
fined close to the hydrostatic core due to the large 
ratio of the impact pressure of the in-falling mat- 
ter with respect to the gas pressure: pv 2 /p = 
7M 2 w 9—12, where pv 2 is the impact pres- 
sure, p is the gas pressure, 7 = 5/3 is the adi- 
abatic coefficient, and M is the Mach number 
of the flow. The density at the shock is low, 
Pshock ~ (1.4 — 5.2) x 10 5 g cm^ 3 , and the temper- 
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ature is high, 



Tqa = T g / (1 + 0.067T 9 ) , 



(6) 



T 



hock 



3 n 
16 kN a 



u 2 « (0.5- 1) x 10 9 K, (3) 



where is the mean molar mass (p = 1.75 g moU 1 
for completely ionized carbon and oxygen matter), 
TVa is Avogadro's number, and k is the Boltzmann 
constant. The accretion shock evolution is there- 
after driven by two opposite effects: inertial con- 
finement due to mechanical energy deposition by 
in-falling matter vs pressure build-up due to nu- 
clear energy release (not included in the present 
simulations). Eventually, the nuclear energy re- 
leased by the burning of a mass M pr i becomes large 
enough for pressure build-up to take over the me- 
chanical energy deposition, causing the accretion 
shock to expand and decouple from the underly- 
ing dense hydrostatic core. If this primer mass is 
large enough to ignite a detonation, then nuclear 
processing of most of the core will be ensured, oth- 
erwise nuclear burning will cease shortly after due 
to adiabatic cooling of the core. 

The rate of mechanical energy deposition at 
the accretion shock, £ mcc h = 27rr 2 /?ov 3 , stays in 
the range (0.5 — 1) X 10 49 erg s _1 during approx- 
imately At ~ 0.3 s. Assuming that nuclear burn- 
ing proceeds just to 28 Si (because p < 10 7 g cm' 3 ) 
the specific nuclear binding energy released is q ~ 
5.8 x 10 17 erg g _1 . An estimate of the maximum 
primer mass that can be accumulated before ex- 
pansion of the accretion shock can be obtained by 
equating the nuclear energy release to the mechan- 
ical energy deposited by in-falling matter M pi iq « 
e mcch Ai, which gives M pri w (2.6 — 5.2) x 10 30 g. 

The radius at which a detonation could ignite 
can be estimated independently from the above 
calculations in the following way. A necessary con- 
dition for a detonation to be initiated is that the 
nuclear timescale must be lower than the hydrody- 
namical timescale, Thydro = 446/^/p s hock ~ 0.6 — 
1.2 s. The nuclear timesc ale of carbon burning. 
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1987tlKhokhlovlll989l) : 



can be calculated as (Blinnikov & Kho khlov 



where is the Frank-Kamcnctskii factor, 



9 = 



3T„ 



(4) 



(5) 



Q = 84.165, and T 9 = T/10 9 K. Finally, r, is 
obtained as: 



c c T 



pq c A T9 Yg exp (q/tI 



/3 
9A 



(7) 



where q c 



4.48 x 10 18 erg mol -1 , c is the 



specific heat, A T9 = 8.54 x 10 26 T 9 5 ^ 6 T 9 3/2 , and 
Yq = 0.5/12 mol g" 1 is the 12 C molar fraction (we 
assume a 50%-50% by mass carbon-oxygen com- 
position). 

Using the above equations with the values 
of Pshock and Tshock from Table [2] the nuclear 
timescale is larger than the hydrodynamical 
timescale, hence a detonation cannot begin just 
behind the accretion shock. However, the condi- 
tion r nuc < Thydro can be met during the journey 
of the shocked matter along its path towards the 
surface of the hydrostatic core. In order to obtain 
an analytical estimate of the distance from the ac- 
cretion shock at which the timescale requirement 
can be fulfilled we need to make some simplifying 
assumptions concerning the shocked flow: 

• Steady state during times small compared 

with Thydro- 

• Spherical symmetry and radial flow (we ig- 
nore any instability in the flow). 

• Optically thick matter (we assume thermal 
equilibrium with radiation). 

• Adiabatic evolution (we ignore the nuclear 
energy input and any energy transfer mech- 
anism, plus optically thick matter which im- 
plies that radiative cooling is inefficient). 

Under these assumptions, the structure of the 
shocked flow in between the accretion shock and 
the core surface can be obtained by solving the 
following set of equations: 
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GM 
r 

constant , 
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constant . 



(8) 

(9) 
(10) 



The starting point for the integration of these 
equations is just the physical state behind the 
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shock (see Table [2]): p sri ock, and T sri0 ck- The 
hydrodynamical and nuclear timescales obtained 
from the solution to the above equations for model 
DF18 are plotted in Fig. [TQ1 where it can be seen 
that carbon burns in less than a hydrodynamical 
time at a radius of ~ 3 000 km. The density and 
temperature at that point are given by the cross 
symbol in Fig. QTJ 

As shown in Table [H the distance from the ac- 
cretion shock to the detonation ignition point is 
sensitive to the mass burned subsonically. The 
Lagrangian mass at which a detonation could be- 
gin, Mdeto, goes from 1.03 M Q for model DF26 
to 1.29 M for model DF11, which is 0.06 M Q 
to 0.03 Mq inwards from the accretion shock, re- 
spectively. The value of Mdcto quoted in Table [2] 
for model DF29, M dcto = 0.38 M Q is far too low 
to be minimally realistic. However, we stress that 
the analytic calculations of the formation of a det- 
onation presented in this section are just approx- 
imations because of the many simplifications in- 
troduced, particularly the assumption of spher- 
ical symmetry and steady state. In the three- 
dimensional picture there are to be expected some 
deviations from the results of these analytic cal- 
culations. 

3.2. Detonation ignition conditions 

Small-scale simulations are needed to determine 
which arc the suitable conditions for detonation 
i gnition in white dwarf matt er. 



Arnett k Livne 



1994ft: iNiemeyer k Wooslevl (|1997ft : iRopke et al 
2007bh performed such studies following the hy- 



drodynamical and nuclear evolution of a uniform 
density microsphere. This thermonuclear bomb 
was lit using a primer consisting in a thermal pro- 
file characterized by a central peak temperature 
and a constant negative thermal gradient down to 
10 8 K at a Lagrangian mass coordinate M pT1 , be- 
yond which the temperature was uniform. Once 
ignited, the sphere was subject to two opposite 
effects: it expanded and cooled due to the ex- 
cess of pressure derived from the sudden inciner- 
ation of the hot primer but, at the same time, 
the released nuclear energy helped keeping a high 
value of temperature. Which one of the two ef- 
fects wins is what determines if a detonation is 
successfully launched or not. These numerical ex- 
periments are therefore characterized by three pa- 
rameters: the central temperature, T c , the uni- 



form sphere density, p, and the mass of the primer, 
M pri . A larger primer mass implies a larger en- 
ergy reservoir and a stronger piston effect, making 
it easier to build a detonation wave. The detona- 
tion initiation conditions i n pure C-0 matter (50% 
each by mass) obtained bvlArnett k Livne (Il994 ) ; 
iNiemeyer k Wooslevl (fl997ft : lRopke et all (|2007bh 
are summaryzed in Table [3] Their results can be 
interpreted either as the minimum primer mass 
for which a detonation is initiated in C-0 mat- 
ter at a given density and peak temperature, or 
as the minimum temperature for which a detona- 
tion is obtained at a given density and mass of the 
primer. In any case, the qualitative trend is that 
the smaller the density the harder is to initiate a 
detonation while, at fixed density, the higher the 
peak temperature the lesser primer mass is needed 
to detonate the sphere. From Table[3]it stems that 
at densities in the range 10 7 — 10 8 g cm -3 a tem- 
perature of ~ 2 x 10 9 K is enough to initiate a 
detonation. At lower densities the primer mass 
necessary to initiate a detonation at 2 x 10 9 K be- 
comes a non- negligible fraction of the star mass, 
hence achievement of detonation initiation condi- 
tions is much more difficult. At p — 3x 10 6 g cm" 3 
and T = 2.3 x 10 9 K a detonation was obtained 
for M pr j = 2 x 10 28 g, implying a thermal gradient 
VT = 200 K cm" 1 , of the same order of the max- 
imum thermal gradient estimated by using Eq. [T] 
VT max = 5 — 150 K cm" 1 . 

The results of the small-scale simulations al- 
low us to look for regions in the three-dimensional 
hydrodynamic simulations that meet some gross 
criteria for detonation ignition. A comparison of 
the results of the present three-dimensional simu- 
lations of white dwarf pulsation with the detona- 
tion ignition conditions summarized in Table [3] is 
shown in Fig. [TT] where we plot in the p-T plane 
the evolution of the fuel particles that achieve 
the most favorable conditions for detonation ig- 
nition, i.e. those particles that attain the max- 
imum temperature after being shocked. Shown 
in this Fig. are the results for five models for 
which M,jefl spans the range from 0.11 M Q to 
0.29 M Q . The two big dots belong to the two 
first rows in Table [3J which are used to delimit 
the p-T conditions for which a C-0 detonation 
is the most probable outcome. As can be seen 
in this Fig., such conditions are clearly reached 
in all the three-dimensional pulsating models ex- 
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cept for DF29, characterized by a mass burned 
subsonically close to the mass needed to unbind 
the white dwarf (0.30 M Q ). The path followed 
by the particles belonging to the different mod- 
els is nearly identical, and very close to adiabatic 
evolution (dot-dashed line in Fig. [TT]) . The point 
at which model DF18 (green points) crosses the 
line defined by r nuc = Th y dro is in good agreement 
with the estimate based on Eqs. IMTOl (dark- green 
cross). In all cases but model DF29, the C-0 deto- 
nation ignition line is crossed at ~ 3 x 10 6 g cm~ 3 . 
The peaks of temperature and density span a 
narrow range: T max w (2.4 — 2.9) x 10 9 K and 
Pmax ~ (4 — 8) x 10 6 g cm~ 3 , high enough to con- 
fidently assume that a CDI will be produced in 
these cases. 

In Fig. [T3 there are shown the thermodynamic 
and chemical structure in a slice of model DF18 
in the neighborhood of the hottest shocked fuel 
particle before (top row) and after (bottom row) 
it reaches detonation ignition conditions (T ~ 
2.3 x 10 9 K, p ~ 3 x 10 6 g cm' 3 ). The deto- 
nation conditions are first reached in a hot spot 
of radius ~ 200 km (mass ~ 10 29 g) affect- 
ing a region with non-uniform chemical compo- 
sition, where the fuel mass fraction is in the range 
X(C + 0) - 0.4- 0.8. A rough estimate of 
the thermal gradient across the hot spot gives 
VT - 10 9 K/200 km = 50 K cm" 1 . Thus, al- 
though the situation is much more complex in our 
models than that addressed in the small-scale sim- 
ulations, we can confidently expect that a stable 
detonation will be obtained a few tenths of a sec- 
ond after the formation of the accretion shock. 
Note that the hot spot that appears close to the 
center in the top row image is due to the presence 
of a clump of ashes, which avoids a detonation to 
be initiated at its position. Note also that, be- 
cause we did not allow for nuclear reactions in the 
present models once the deflagration phase ended, 
the chemical changes between both snapshots are 
only due to advection of particles of different com- 
position. 

The maximum temperature and density reached 
in model DF29, T max = 1.42 x 10 9 K and p max = 
1.3 x 10 6 g cm -3 , are seemingly too low fo r igniting 
a stable detonation. iRopke et alJ (l2007bh did not 
obtain a stable detonation for p = 10 6 g cm -3 and 
a temperature of T = 3.0 x 10 9 K even for a primer 
mass as high as M pr j = 3.0 x 10 30 g. This mass is 



comparable to the maximum primer mass derived 
in Sect. O for model DF29: M pri = 2.6 x 10 30 g 
(see also Table HJ). 

For comparison purposes we show in Fig. [11] 
the path followed during the recontraction phase 
by the PDD model described in Section 12.21 but 
with the nuclear reactions turned off. The dis- 
played path is that drawn by a fuel particle in the 
neighborhood of the quenched flame. The path 
of the particle penetrates the C-0 detonation re- 
gion at p ~ 7 x 10 7 g cm" 3 and T ~ 1.2 x 10 9 K, 
and afterwards reaches a maximum temperature of 
T max = 1.5 x 10 9 K. However, note that this region 
of the p-T plane has not been thoroughly explored 
in order to determine the conditions suitable for 
detonation ignition (see Table [3] and references 
therein). Thus, in this region the threshold line 
is just an extrapolation from lower density points. 
It is therefore clear that the conditions achieved 
during the re-contracting phase of PDD models 
are quite different from those obtained with PRD 
models. 

3.3. Resolution issues 

Our numerical resolution is much coarser than 
the expected detonation width at the densities 
of detonation ignition. Thus, we have explored 
the sensitivity of the maximum temperature and 
density reached during white dwarf pulsation to 
the resolution of the numerical models. To this 
end we have computed additional simulations of 
the pulsating phase of models DF18 and DF29 
with increased resolution. The starting point of 
these additional models was the time of maxi- 
mum expansion during the pulsation. The par- 
ticles that fulfilled certain selection criteria were 
splitted into N 8p children particles, and the star 
evolution was again followed with the same three- 
dimensional SPH code described before (the term 
particle splitting was used fo r the first time by 
Kitsionas fc Whitwortb l2002h . In order to save 
computer time, splitting was confined to the neigh- 
bors of the shocked fuel particles that reached the 
maximum temperature (Fig.lTTjl. and to the neigh- 
bors of their neighbors. We performed simulations 
for N sp =2, 4, 10, and 100. In every case the origi- 
nal particle mass was shared evenly between their 
children, which resulted in a maximum mass res- 
olution of 1.1 x 10 26 g. The initial position of the 
children particles was randomly selected within a 
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sphere centered on the position of the parent parti- 
cle and radius given by its smoothing length. They 
were therefore assigned the same temperature and 
chemical composition as the parent particle, while 
their velocity was obtained by interpolation of the 
velocity field at the position of the children. 

The results of our resolution study are shown 
in Table 01 where there are given the maximum 
density and temperature reached by children par- 
ticles as well as their maximum spatial resolution, 
Aa; m i n . Even though the maximum temperature 
and density for a given M^cH show some scatter 
with increasing resolution they all stay at the same 
side of the C-0 detonation line. The departure 
from the conditions obtained with our standard 
non-splitted models is quite modest and does not 
affect our conclusions: all the variants of model 
DF18 fulfilled the conditions for a C-0 detonation, 
while all the variants of model DF29 did not. 

In our highest resolution simulations of model 
DF18, the increase in resolution does not lead to 
a too steep thermal gradient that might compro- 
mise the achievement of detonation conditions. In 
model DF18spl00, the primer mass (calculated 
as the total mass of fuel particles with p > 3 x 
10 6 g cm" 3 and VT < 200 K cm" 1 that are neigh- 
bors of the hottest children) rises to 1.8 x 10 28 g. 
This mass is similar to the minimum mass required 
to detonate C-O, as given in Table [3] for the same 
density but for a peak temperature substantially 
lower than the maximum temperature achieved in 
model DF18spl00. As we have not been able to 
achieve convergence of the maximum temperature 
and density, it is not clear if further increase of the 
resolution might result in a thermal gradient too 
steep to launch a detonation. 

4. A helium cap on the white dwarf? 

One way to favor a detonation ignition in C-0 
matter is to mix it with a small amount of He, 
which both is more reactive and releases more nu- 
clear energy than carbon. To allow such a mix- 
ing, some mass of He must be present in the 
white dwarf at the moment of thermonuclear run- 
away. Indeed, a very small cap made of He might 
rest atop of the white dwarf as a result of pre- 
vious accretion from the companion star in the 
binary system. This helium could either have 
been accreted directly from a He star or result 



from the nuclear processing of accreted hydro 



gen. Hoflich fc Khokhlov ( 19961 ) suggested that 



the presence of a ~ 0.01 M Q cap of He would im- 
prove the match between calculated and observed 
light curves. This amount of He might come from 
the surrounding accretion disk if the accretion rate 
were not constant, for instance if it dropped at 
som e time below the critical rate for steady burn- 
ing. Cumming et alJ (|l994bl |al) detected lines coin- 
cident with two He I lines at 2.04 /mi and 1.052 /im 
in the spectra of the spectroscopically normal SN 
1994D. The absorption feature at 1.05 /xm has 
been since dete cted in the early spectra of many 
other SNIa (see Nomoto et al.ll2003l Pi gnat a et al 



2008), although its identification is problematic 



because that line could rather be due t o Mg II 
(jWheeler et alJll998t lHatano et aT1ll999l) . How- 
ever, as pointed out by iMazzali fc Lucy ( 1998 ). 
He I lines may be blended with magnesium lines 
as well as with lines of other intermediate-mass 
elements. Hence, up to ~ 0.01 M of He might 
remain undetected in SNIa, hidden by stronger 
magnesium lines. 

If there was actually such small cap of he- 
lium, it would have several chances of mixing with 
the C-0 outer layers of the white dwarf. The 
first opportunity would arise during the convec- 
tive simmeri ng phase shortly before the hydrody- 
namic event dKuhlen et alj|2006t iPiro fe Bildstenl 
20081: IChamulak et all l2008l )~Even though con- 
vection during simmering is not expected to reach 
the white dwarf surface, it can excite pulsating 
and non-spherical mode instabilities of the star 
that might allow diffusion of He into the un- 
derlying C-0 layers. A second chance for mix- 
ing He would take place during the deflagra- 
tion phase of the explosion itself. Turbulence 
and, especially, Rayleigh- Taylor instability induce 
large radial excursions of huge volumes of ashes 
that destabilize the whole white dwarf, and can 
drive the mixing of an eventual external He layer 
with the C-0 beneath it. Hence, even though it 
may seem a speculative scenario, one cannot dis- 
card the possibility that He is present in some 
small mass fraction within the matter that is 
prone to detonate after a white dwarf pulsation. 
Mixing of He within a detonating volume seems 
more likely in scenarios in which the transition 
to detonation develops in the outermost layers 
of the white dwarf, e.g. those CDI described in 



10 



Plewa et all (12004 . iBravo fc Garcia-Senj (|2006f) 
and lRopke et all (|2007bf )~ 

We have explored the conditions for stable deto- 
nation ignition in white dwarf matter composed by 
C-0 contaminated with He, following a method- 
olog y similar to that of Arnett fc Livnel ( 1994 ) 
and iNiemever fc Wooslevl (|1997l h Using a one- 
dimensional code we have followed the hydrody- 
namical and nuclear evolution of a uniform den- 
sity sphere made of equal amounts of carbon 
and oxygen with a small fraction of He. The 
setup of these numerical experiments, as well as 
the physi cs included, were similar t o those de- 
scribed in lNiemeyer fc Wooslev ( 1997b . These nu- 
merical experiments are characterized by four pa- 
rameters: the central temperature, T c , the uni- 
form sphere density, p, the mass of the primer, 
M pr i, and the He mass fraction, JT(He). In or- 
der to limit the dimensionality of the problem 
we have performed calculations at a fixed den- 
sity: p — 1.5 x 10 6 g cm~ 3 , while we have ex- 
plored the following ranges of the rest of parame- 
ters: 1.0 x 10 9 K T c < 1.8 x 10 9 K, 3 x 10 25 g 
M pri ^ 3 x 10 29 g, 0.02 X(Re) ^ 0.20. Be- 
sides, two different He distributions have been ex- 
plored: either it was concentrated in the primer 
volume, or it was uniformly distributed through- 
out the whole microsphere. For the extreme values 
of primer mass considered the ball sound crossing 
time ranges from 0.01 s (M pri = 3 x 10 25 g) to 0.3 s 
(Af pri = 3 x 10 29 g). 

Our results, summarized in Tabled show that 
addition of a small quantity of He makes deto- 
nation initiation far easier than in pure C-O. In 
the numerical experiments in which He was con- 
strained to the primer its role was limited to pro- 
vide a larger nuclear energy release as compared to 
pure C-O, while distributing it uniformly through- 
out the sphere allowed also to take advantage of 
the larger reactivity of He. Table [5] shows that 
with a He mass fraction of 10% a detonation is 
feasible even at a peak temperature as small as 
10 9 K provided that a large enough primer mass 
(> 1.5 x 10~ 4 M Q ) is present. Note that in 
this case the total mass of He needed to deto- 
nate the microsphere is just > 1.5 x 10~ 5 M Q . 
For a 5% He mass fraction the minimum tempera- 
ture needed to achieve a stable detonation rises to 
1.8 x 10 9 K for the same primer mass. The calcu- 
lations with a peak temperature of 1.4 x 10 9 K and 



M pr i = 3x 10 g show interesting results: a 20% 
mass fraction of He constrained to the primer was 
not able to initiate a detonation, while a modest 
10% He mass fraction uniformly distributed was 
enough to detonate the microsphere. These re- 
sults suggest that the main reason by which Ho 
facilitates detonation ignition in C-O matter is be- 
cause of its large reactivity rather than because of 
its large nuclear energy release. 

The application of the present results to the 
question of the feasibility of a CDI in model DF29 
is inconclusive. Even disregarding the speculative 
nature of the presence of a He cap onto the white 
dwarf, the small Lagrangian mass at which the ac- 
cretion shock forms in model DF29, ~ 0.38 M Q , 
poses another difficulty to mixing He down to 
such a deep place. In spite of these uncertain- 
ties, we have simulated the detonation phase of 
model DF29 contaminated with He, whose re- 
sults will be presented in detail in paper II: for 
MdeB — 0.29 M & we have not been able to obtain 
a stable detonation irrespectively of the presence 
of any reasonable amount of He on top of the white 
dwarf. 

5. Conclusions 

The initiation of a detonation is a highly non- 
linear process in which a modest variation in the 
environmental conditions can have a large im- 
pact on the outcome. Theoretical analysis, small- 
scale simulations and large-scale numerical sim- 
ulations of white dwarf explosions powered by 
slow subsonic flames do not favor presently a 
deflagration-to-detonation transition during the 
expansion phase of the explosion, although it is 
not discarded either. The exception may be DDT 
induction by transition of th e flame to the dis- 
tributed combustion regime (jRopke fc Niemever 
20071 : 1 Wooslevl 120071 ). but there remain many un- 
certainties about its feasibility in an explodin g 
white dwarf (jLisewski et al.|[2000l|Pan et al.ll2008l ). 

We have shown that a confined detonation igni- 
tion can be induced by an accretion shock, formed 
around a carbon-oxygen core after WD pulsation 
following subsonic burning of a mass Md c fi J$ 
0.3 M Q . The mechanism of CDI is robust, and 
detonation ignition conditions are reached non- 
marginally within the shocked flow for a wide 
range of initial configurations. However, if the nu- 
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clear energy released before the pulsation is close 
to the white dwarf binding energy the accretion 
shock is not efficient enough to detonate the mix- 
ture of carbon and oxygen. 

We have also addressed the effects that the pres- 
ence of a He cap on top of the white dwarf, at the 
time of thermal runaway, could have for detona- 
tion ignition. Such He cap may partially mix with 
the underlying carbon-oxygen layers during differ- 
ent phases of the evolution of the white dwarf. 
Through small-scale one-dimensional spherically- 
symmetric calculations of detonation ignition in 
carbon-oxygen matter contaminated with traces 
of He, we have determined the conditions under 
which the He presence can allow the formation of 
a detonation. It turns out that a local contami- 
nation by a 10% mass fraction of He is enough to 
allow detonation ignition for a peak temperature 
of just 10 9 K, provided that the thermal gradi- 
ent is shallow enough (implying a large enough 
primer mass). If the He mass fraction is reduced 
to A (He) < 0.05 the peak temperature necessary 
to detonate rises to <~ 1.8 x 10 9 K. These results 
might be relevant not only to CDI in PRD models, 
but also to the more general problem of detonation 
ignition in any delayed detonation scenario. 

Our results show that the PRD scenario of 
white dwarf delayed detonation is feasible and 
worth pursuing its study. PRD models may ex- 
plain many observational features of a wide range 
of SNIa events, but probably they cannot cover 
the whole range of energy and radioactive masses 
deduced from observations. What we have shown 
in the present work is that the PRD scenario pro- 
vides a realistic way by which white dwarfs can 
explode. In paper II we will analyze the explosion 
phase and compare the final outcome with obser- 
vational features of typical SNIa. 
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Table 1 

Initial size and position of the hot bubbles' 1 







z c 


M b 


(km) 


(km) 


(km) 


(g) 


-56.5 


74.5 


87.0 


1.59 x 10 30 


101.3 


-87.7 


-52.2 


1.48 x 10 30 


106.5 


41.3 


47.2 


1.53 x 10 30 


115.9 


180.3 


-195.0 


1.09 x 10 30 


-145.8 


-83.4 


35.8 


1.44 x 10 30 


275.9 


-101.5 


-80.4 


0.87 x 10 30 


32.4 


239.6 


-61.6 


1.35 x 10 30 



a In this table x c , y c and z c are the coor- 
dinates of the center of the bubbles, while 
Mb gives the mass of each bubble 
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Table 2 

Physical conditions of shocked matter 



Model 


DF11 


DF18 


DF29 


Mdefl (M ) 


0.11 


0.18 


0.29 


tacc (s)" 


4.00 


7.18 


12.17 




12 


12 


9 


Pshock (g cm~ 3 ) 


5.2 x 10 5 


3.3 x 10 5 


1.4 x 10 5 


T 3hock (K) 


8.1 x 10 s 


10° 


4.5 x 10 8 


Smeoh (erg s _1 ) 


1.0 x 10 49 


10 49 


5 x 10 48 


M pri (g) 


5.2 x 10 30 


5.2 x 10 30 


2.6 x 10 30 


Thydro (s) b 


0.62 


0.78 


1.21 


Mdeto (M )° 


1.29 


1.10 


0.38 


tdoto (s) d 


4.3 


7.5 





a Time at which the accretion shock is formed 

b Hydrodynamic timescale just behind the accretion shock 

c Lagrangian mass at which detonation ignition conditions 
would be achieved. See Sect, [3TT1 

d Time at which detonation ignition conditions arc achieved 
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Table 3 

Summary of detonation initiation conditions in C-O matter 3 - 



p 




Afpri 


(g cm" 3 ) 


(10 9 K) 


(g) 


3 x 10 6 


2.3 


2 x 10 28 


10 7 


1.9 


1.5 x 10 27 


10 7 


2 


1.1 x 10 24 


in 7 


2.2 


2 x 10 25 


10 7 


2.8 


2.5 x 10 23 


10 7 


3.5 


2 x 10 23 


2 x 10 7 


2 


2.9 x 10 21 


3 x 10 7 


2 


4.9 x 10 20 


3 x 10 7 


5.2 


2 x 10 19 


5 x 10 7 


2 


1.1 x 10 19 


10 s 


2 


5.5 x 10 17 


10 s 


6.2 


2 x 10 15 



a Compiled from |Arnett &: Livnel 
|[7994T); iNiemever Wooslevl H1997I) ; 
IRopke et al.l I l2u07bl) 
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Table 4 

Maximum temperature reached with different mass resolutions 



Model 


Pmax 

(g cm- 3 ) 


Tmax 

(10 9 K) 


Alllmin* 
(g) 


(km) 


DF18 


4.6 X 10 6 


2.40 


1.1 x 10 2S 


83 


DF18sp2 


4.5 X 10 6 


2.36 


5.5 x 10 27 


66 


DF18sp4 


2.7 x 10 6 


2.76 


2.8 x 10 27 


63 


DF18splO 


4.6 x 10 6 


2.74 


1.1 x 10 27 


39 


DF18spl00 


3.7 x 10 6 


3.70 


1.1 x 10 26 


19 


DF29 


1.3 x 10 6 


1.42 


1.1 x 10 2S 


126 


DF29splO 


2.6 x 10 6 


1.91 


1.1 x 10 27 


47 


DF29spl00 


1.7 x 10 6 


1.95 


1.1 x 10 26 


25 



a Maximum mass resolution 

b Maximum spatial resolution, given by Ax-_ 
(3Am min /47rp max ) 1/3 
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Table 5 

Detonation initiation in C-O matter with traces of He 



T 


M pri a 


X(Ho) b 


Hc-distribution c 


Detonation? 


VT 


(10 9 K) 


(g) 








(K cm" 1 ) 


1.0 


3 x 10 29 


0.10 


primer 


y 


28 


1.4 


3 x 10 25 


0.10 


unif. 


y 


830 


1.4 


3 x 10 27 


0.02 


primer 


n 


180 


1.4 


3 x 10 27 


0.05 


primer 


n 


180 


1.4 


3 x 10 27 


0.10 


primer 


n 


180 


1.4 


3 x 10 27 


0.20 


primer 


n 


180 


1.4 


3 x 10 27 


0.10 


unif. 


y 


180 


1.4 


3 x 10 29 


0.05 


primer 


n 


39 


1.4 


3 x 10 29 


0.10 


primer 


y 


39 


1.8 


3 x 10 27 


0.05 


unif. 


n 


230 


1.8 


3 x 10 29 


0.05 


primer 


y 


50 



a In all calculations the density was p — 1.5 X 10 6 g cm -3 
b He mass fraction 

c Primer: He concentrated in the primer. Unif.: He uniformly distributed 
through the whole microsphere 
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0.1 1 10 

time (s) 

Fig. 1. — Evolution of the models in Table 
through the deflagration and pulsation phases. 
Shown are the evolution of the burnt mass, Mb 
(left axis scale), and the central density (right 
axis scale), for models DF11 (dashed lines), DF18 
(solid lines), and DF29 (dotted lines). Once the 
white dwarf reached its minimum central density 
the thermonuclear burning virtually had quenched 
(t > 2 s), and the hydrodynamic evolution was 
thereafter followed with the same SPH code but 
with the nuclear reactions and the deflagration 
propagation algorithm switched off. The two ver- 
tical dash-dotted lines mark the times of the snap- 
shots shown in Fig. [12] (model DF18) 




Fig. 2. — Evolution of model DF18 during the first 
second of the hydrodynamical simulation (defla- 
gration phase). Shown are the isosurfaces defined 
by an ash mass fraction X as h = 0.5, and volume 
renderings of the density (in blue shades). The 
minimum density depicted is 10 6 g cm -3 in the 
two first snapshots, 3 x 10 5 g cm~ 3 in the third 
snapshot, and 3 x 10 4 g cm -3 in the last one. The 
coordinate origin is at the center of mass of the 
white dwarf. The axes units are 1 000 km 
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1 2 3 4 5 6 7 

time (s) 



Fig. 3. — Evolution of the white dwarf during the 
first pulsation of model DF18. Each curve shows 
the radius of a sphere that encloses a constant 
mass. Starting from the center outwards the mass 
of each curve increases by 0.0345 M Q up to the 
thick curve (which belongs to ta Lagrangian mass 
of ~ 1.10 M Q ), from which the increase is reduced 
to 0.0138 M . The accretion shock forms approx- 
imately 4-5 shells above the thick line. The two 
vertical dash-dotted lines mark the times of the 
snapshots shown in Fig. [12] (model DF18) 




Fig. 4. — Map of the fuel (C+O) mass fraction 
in and around the dense core during the first pul- 
sation of model DF18. The three surfaces shown 
in the pictures belong to the planes x = 0, y = 
5 000 km, and z = 5 000 km. The axes units are 
1 000 km 




Fig. 5. — Thermal and mechanic structure of the 
white dwarf at the time of accretion shock forma- 
tion in model DF18 (t ~ 7.2 s). The three sur- 
faces shown and the axes units are as in Fig. 0] 
Here, as in Figs. [5] and 2J each three-dimensional 
image represents a scalar field reconstructed from 
the computed SPH particles properties (mass, po- 
sition, velocity, smoothing length, chemical com- 
position, temperature) through standard SPH in- 
terpolation into a regular Cartesian lattice 
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Fig. 6. — Distribution of specific energy within 
the white dwarf during the first pulsation in a 
one-dimensional calculation of a pulsating delayed 
detonation model (model PDDe in Badenes et al. 
2003) vs the three-dimensional model DF18: spe- 
cific thermal energy (dotted-line) and kinetic en- 
ergy (dashed- line) . The top panel shows profiles 
obtained from angle- averaged versions of the DF18 
model at a time t=4.2 s after initial runaway. The 
bottom panel shows the profiles of model PDDe 
at a similar elapsed time 
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Fig. 7. — Same as previous Fig., but for the chem- 
ical composition. Solid-line: C+O, dotted-line: 
products of the incineration to nuclear statistical 
equilibrium (mainly Fe-Ni) 




0.5 1 
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Fig. 8. — Evolution of velocity (solid lines, in 1,000 
km s -1 ) and entropy (dotted lines, in fee per par- 
ticle) calculated with a one-dimensional code af- 
ter an angle-averaged version of the DF18 model 
starting at t = 4.2 s. Times plotted arc 4.2 s, 
6.2 s, 7.0 s, and 8.2 s. For clarity, neither the den- 
sity profile at 4.2 s nor the entropy profiles at 4.2 s 
and 6.2 s have been labelled 
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Fig. 9. — Evolution of velocity of a PDD model 
that burns subsonically the same mass as model 
DF18. Times plotted are 4.2 s, 6.2 s, 7.0 s, 8.2 s, 
and 11.4 s. The dotted line gives the escape ve- 
locity at t = 4.2 s 
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Fig. 11. — Maximum temperature and density 
reached during the pulsating phase of several mod- 
els. The colored points show the track of the 
hottest shocked fuel particles of each model. The 
shaded area is the region where a C-0 detonation 
can be successfully launched, extrapolated from 
the t wo first rows in Table [3] ( Ropke fc Niemever 
2007t ). The dash-dotted line represents a T oc p 1 / 3 
relation, plotted for reference. The dotted line 
shows the location of r nuc = Th y dro- The point at 
which this line is crossed in the DF18 model ac- 
cording to the solution of Eqs.[8TfTUlis given by the 
dark-green cross symbol. The time at which differ- 
ent models achieve detonation ignition conditions 
(at the entrance to the shaded area) is given in the 
last row of Table [5] The SPH models summarized 
in this figure did not include nuclear reactions 



Fig. 10. — Comparison of the hydrodynamic 
timescale (solid line) and the timescale of car- 
bon burning (dashed-line) within the shocked flow. 
The zero of the horizontal axis marks the location 
of the accretion shock 
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Fig. 12. — Thermodynamic structure (left) and 
fuel mass fraction (right) in a slice of model DF18 
in the neighborhood of the hottest shocked fuel 
particle plotted in Fig. [TlJ The top row belongs 
to a time shortly before the formation of the accre- 
tion shock, while the bottom row belongs to the 
time at which the particle reaches C-0 detonation 
conditions (shaded area in Fig. [TT|) . i.e. ~ 0.3 s 
later than the top row (see Table [H and Figs. Q] 
and [3]) . The position of the particle is shown by 
a black dot located close to the bottom and to 
the center of each slice. The temperature is given 
in units of 10 9 K, while the density contours be- 
long to, from the center outwards and in units of 
g cm -3 , logp = 8, 7.5, 7, 6.5, 6, and 5.5 (the con- 
tour belonging to logp = 6.5 is labelled in the bot- 
tom row slices for clarity). The images have been 
obtained from the computed properties of the SPH 
particles as explained in Fig. [5j The axes units are 
1 000 km and their directions are the same as in 
Figs. H 12 and [5] 



